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ABSTRACT 

HH 212 is a nearby (400 pc) highly collimated protostellar jet powered by a Class 0 source in Orion. We 
have mapped the inner 80" (~ 0.16 pc) of the jet in SiO (J = 8 — 7) and CO (J = 3 — 2) simultaneously 
at ^ 0"5 resolution with the Atacama Millimeter/Submillimeter Array at unprecedented sensitivity. 
The jet consists of a chain of knots, bow shocks, and sinuous structures in between. As compared to 
that seen in our previous observations with the Submillimeter Array, it appears to be more continuous, 
especially in the northern part. Some of the knots are now seen associated with small bow shocks, 
with their bow wings curving back to the jet axis, as seen in pulsed jet simulations. Two of them are 
reasonably resolved, showing kinematics consistent with sideways ejection, ossibly tracing the internal 
working surfaces formed by a temporal variation in the jet velocity. In addition, nested shells are seen 
in CO around the jet axis connecting to the knots and bow shocks, driven by them. The proper motion 
of the jet is estimated to be ^ 115T50 km s -1 , comparing to our previous observations. The jet has a 
small semi-periodical wiggle, with a period of ~ 93 yrs. The amplitude of the wiggle first increases with 
the distance from the central source and then stays roughly constant. One possible origin of the wiggle 
could be the kink instability in a magnetized jet. 

Subject headings: stars: formation — ISM: individual: HH 212-ISM: jets and outflows. 


1. INTRODUCTION 

Protostellar jets represent one of the most intrigu¬ 
ing signpo sts of star formation (for recent review see 
iFrank et al.l 120141 ). They are highly supersonic and col¬ 
limated, and some with wiggles in their trajectories. They 
are seen with knotty and bow-like shock structures. They 
can be launched from accretion disks around protostars, 
allowing us to probe the accretion process closest to the 
protostars. Two c ompetin g yet similar MHD models, 
the X-wind model (IShu et al.l 120001 and disk-wind model 
([Konigl fe Pudrit3 l2000lh have been proposed to launch 
these jets. Preliminary measurements of specific angular 
momentum of some jets are consistent with thes e two mod¬ 
els (see, e.g. dLee et al.ll2QQ8[ ICoffev et al.ll201ll) . Polariza¬ 
tion obs ervations also suggest that the j ets could be mag - 
netized ([Carrasco-Gonzalez et al.l[20101: iLee et al.ll2014bl h 
as predicted in the two models. Therefore, more detailed 
studies of the jet properties are needed to determine the 
jet launching model, and thus the accretion process closest 
to the protostars. 

Molecular outflows are shell-like structures surrounding 
the jets and thus can be used to probe the physical prop¬ 
erties of the jets indirectly. They are believed to be out¬ 
flow cavity walls consisting mainly of the ambient material 
swep t-up by the bow shocks in the jets (see, e.g. dLee et al.l 
I2001IL Since some molecular outflows show wide-opening 
cavity walls near the protostars, especially in the late r 
stage of star formation (|Lee fc Holl2QQ5l: lArce et al.ll2013lh 
an additional unseen wide-angle tenuous winds might be 
required. Indeed, current theoretical jet models also pre- 
dict such winds surrounding the highly collimated jets 
(IShu et al.l 120001: iKonigl fc Pudrit j l2000fi . However, more 


observations are still needed to confirm this possibility, es¬ 
pecially in the early phase of star formation. 

This paper is a follow-up study to our previous Submil- 
limeter Arra y (SMA) study of the HH 212 jet and out¬ 
flow system (|Lee et al.ll2QQ7l [20081. presenting our study 
in SiO and CO with Atacama Millimeter/Submillimeter 
Array (ALMA) at an angular resolution of ~ 0"5 with 
more than 10 times higher sensitivity. HH 212 system is 
nearby in the L1630 cloud of Orion at a distance of ~ 400 
pc. Its jet was discover ed in shock excite d H 2 emission 
(IZinnecker, McCaughrean, fc Ravnerl11998 ). It is powered 
by a low-luminosity (~ 9 L©) Class 0 protostar IRAS 
05413-0104. It interacts with the ambient materia l, driv¬ 
ing a collimated CO outflow around it ([Lee et al.l l2000h . 
A rotationally supported disk is believed to have formed 
around the protostar in o rder to launch the jet (ILee et al.l 
12014 ICodella et al.ll2014lh Lying close to the pla ne of the 
sky (< 5°, IClaussen et al.l[19981: ILee et al.ll20Q7h . the jet 
is one of the best candidates to investigate the jet prop¬ 
erties. Here, we refine the proper motion and thus the 
mass-loss rate of the jet. We study the shock structures 
in the jet, the morphological relationship between the CO 
outflow and the jet, and discuss if there is a need of a 
wide-angle tenuous wind around the jet. The jet has a 
small semi-periodical wiggle. We model the wiggle and 
discuss possible origins for the wiggle. Note that the small 
jet rotation tentative ly found at higher resolution of ^ 
0"35 in SiO ([Lee et alJl2008 ) can not be confirmed here at 
a lower resolution of ~ 0"5. 


2. OBSERVATIONS 
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Observations of the HH 212 protostellar system were 
carried out with ALMA on 2012 December 1 during the 
Early Science Cycle 0 phase. The d etai ls of_these obser¬ 
vations have been reported in iLee et al.l ( 2014 ). and thus 
only important information is reported here. A 9-pointing 
mosaic was used to observe the jet in this system within 
~ 40" from the central source. Both CO J=3-2 and SiO 
J=8-7 lines were observed at a velocity resolution of ^ 
0.21 km s -1 per channel. With super-uniform weighting, 
the synthesized beam has a size of 0753x0747 at a po¬ 
sition angle (P.A.) of ~ 20° in SiO, and 0756x0747 at 
a P.A. of 43° in CO. The rms noise level is ~ 7.7 mJy 
beam -1 (i.e., 0.32 K) for the SiO channel maps, and ^ 
7.2 mJy beam -1 (i.e., 0.28 K) for the CO channel maps. 
The velocities in the channel maps are LSR. The systemic 
velocity in this region is assumed to be V sys = 1.7±0.1 km 
s -1 LSR, as in lLee et al. ( 200 7). The central source is as¬ 
sumed to be at the peak position of the continuum, which 
is amoo) = 05 h 43 m 51?4107, J (200 o) = -01°02'53"167, as 
in lT.ee et, a.U (12014 . 

3. RESULTS 

In order to provide readers with a more complete pic¬ 
ture of the jet and outflow system, our ALMA maps 
are plotted together with an earlier H 2 map of the jet. 
Also, for a better view, all the maps are rotated by 22.5° 
clockwise so that the jet axis is aligned with the y-axis. 
The H 2 map shows that the jet consists of a chain of 
knots and bow shocks, and sinuous (continuous) structures 
in between. It was made by combining the 3 observa- 
tionsjn October 2000, October 2001, and January 2002 
([McCaughrean et “all 12002 ), and is thus roughly 11 years 
earlier than our ALMA observations. Since the proper mo¬ 
tion of the jet is estimated to be ~ 0706 per year later here 
in our paper, the H 2 jet must have moved down the jet axis 
by ^ 077, as compared to the jet seen in our ALMA maps. 
Since this position shift is smaller than the size of a knot 
(> 1"), the H 2 map can still be used to show the rough 
positions of the knots, bow shocks, and sinuous structures 
to be compared with our ALMA maps. 

3.1. Nested Shells in CO 

With ALMA, CO emission can now be mapped at higher 
angular resolution and sensitivity. Figure |T] shows the CO 
maps from blueshifted to redshifted velocities (with a ve¬ 
locity interval of 4.23 km s -1 ) in comparison to the H 2 
map. The blueshifted emission is mainly to the north 
and redshifted mainly to the south of the central source. 
Nested shell structures are seen extending to the north 
and south from the central source. As shown in Figure [2j 
the nested shell structures can be better seen in the CO 
maps integrated over selected velocity ranges that show 
the shell structures. In the south (Figure [2^), a wide- 
opening shell is seen opening to the south from the central 
source, connecting to the big bow shock SB3 further out. 
Likewise in the north (Figure^), a similar wide-opening 
shell is seen opening to the north, connecting to the big 
bow shock NB3 further out. These wide-opening shells 
ar e also seen in the lower transition line of CO at J=2- 
1 (|Lee et al.ll2006lf ). likely tracing the outflow cavity walls 
consisting mainly of the ambient material swept up by the 
big bow shocks. At the base around the central source 


(see the zoom-in in Figure®), the shells fit right into the 
bay of the fl atte ned envelope detected in 350 GHz con¬ 
tinuum (|Lee et al.ll2014l ). This morphological relationship 
suggests that the shells, as they expand, are sweeping up 
and thus excavating the flattened envelope at the base. 

Inside these wide-opening shells, internal shells are seen 
connecting to bow shocks NB1/2 and SB 1/2 from the cen¬ 
tral source (Figure [2]). Here internal shel ls are refe rred as 
the shells located inside the cavity walls ([Lee et al.ll2001 ). 
In the north, the internal shell is also affected by the wings 
of bow shocks NK7 and NK8 inside the shell (Figure QJ-g). 
Inside these shells, CO emission is also seen surrounding 
other smaller knots closer to the source, forming narrow 
internal shells around them. These inner shells are bet¬ 
ter seen as we zoom in to the inner part at low velocities 
as shown in Figure [3j For the prominent pair of knots, 
NK1 and SKI, where the shocks are strong, hollow in¬ 
ternal shells are clearly seen connecting^ to them (Figs [3)3 
and d), as seen earlier in ILee et al.l ( 2007 ). In summary, 
nested shells are seen in CO connecting to the bow shocks 
and knots, driven by them. 

3.2. Molecular Jet in SiO and CO 

Figured] shows the maps of the SiO emission and high- 
velocity (HV) CO emission in comparison to th e H 2 map. 
As compare d to the previous SMA maps in ILee et ahl 
(|2007l l2QQ8h . the maps here are created using slightly 
larger velocity ranges, because ALMA detected the emis¬ 
sion in slightly larger velocity ranges at higher sensitivity. 
For the SiO map, we use the velocity from —23.03 to 17.53 
km s -1 , instead of —21.1 to 16 km s -1 . For the CO map, 
low-velocity emission is excluded to avoid the shell con¬ 
tamination. On the blueshifted side, we use the velocity 
from —21.65 to —8.95 km s -1 , instead of —18.4 to —7.1 km 
s . On the redshifted side, we use 9.56 to 18.77 km s -1 , 
instead of 3.4 to 13.3 km s -1 . Notice that the high-velocity 
ranges are not symmetric with respect to the systemic ve¬ 
locity because the jet has a higher mean velocity and a 
larger velocity range on the bl ueshifte d side than the red¬ 
shifted side (see Figure 5 in iLee et all-2007 ), probably due 
to different inclination angles of the jet in the northern 
and southern parts Lee2007. As found in ILee et all ( 2007 . 
l2008l h both SiO emission and HV CO emissions trace the 
jet well, arising from the sinuous structures and knots in 
the jet seen in H 2 . Notice that since the H 2 map was made 
earlier, the SiO and CO emission peaks are seen slightly 
ahead of the H 2 emission peaks. As discussed in lLee et al.l 
( 2007 ). the jet has a slight bending of ^ 1.5° to the west 
and a small semi-periodical wiggle in the trajectory. 

Due to higher sensitivity of ALMA, SiO emission be¬ 
comes better detected toward the cent r al so urce position, 
as compared to that seen in ILee et all (200 8), likely aris¬ 
ing from the jet base near the source as seen injthe lower 
transition line of SiO at J=5-4 ([Codella et alJ[2007h . Also, 
the jet appears to be more continuous in CO, especially 
in the northern part. In addition, three more knots are 
detected along the jet axis in the south, where no clear 
H 2 knot was detected. One is detected in SiO and CO in 
between knots SK5 and SK7, and is thus labeled as knot 
SK6 as a counterpart of knot NK6 in the north. Two are 
detected in CO in between knots SK6 and SK7, linking 
the two knots. The H 2 and CO emission seen in between 
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NK6 and NK7 could be their counterparts in the north. 
For the prominent knot SKI where the shock is strong, 
the SiO emission is seen only at the bow tip, while the CO 
emission is seen mainly in the bow shock wings. The limb- 
brightened bow wings structure of the CO emission can be 
seen at lower velocity ( Figure |3]i). S imila r morphology is 
also seen in L 1157 in iGueth et al.l (1 19981) . probably be¬ 
cause CO is destroyed or excited to higher transition lines 
at the tip where the shock is stronger. 

The jet can also be seen at lower velocity (see Figure 
[Tj) because of the low inclination angle of the jet to the 
plane of the sky. Interestingly, at higher angular reso¬ 
lution and sensitivity, the CO emission at 5.86 km s -1 
appears to show a helical structure along the jet axis in 
the south (see Figure [3ji ) that is not seen before. This 
helical structure is actually consisted of a chain of small 


1.2 Jy beam -1 km s -1 . As discussed in lLee et~ahl ( 2007 ), 
the mean excitation temperature of the CO emission can 
be assumed to be ^ 50 K. Since the brightness tempera¬ 
ture in the jet is mostly < 20 K, the CO emission can be 
assumed to be optically thin. Assuming LTE, the column 
density of CO is estimated to be ^ 2.3 x 10 16 cm -2 . The 
mean jet (H 2 volume) density can be derived from the CO 
column density using the conversion equation, e.g., Eq. 17 
in iLee et al.l (|2014f ). With a CO abundance of 1 8.5 x 10 -5 
and a jet diameter of ~ 0"2 (|Cabrit et al .11 200 71 : ILee et all 
2008), the mean jet density is estimated to be ^ 5.6 x 10 5 
cm -3 , similar to the mean value derived from HCO + in 
ILee et all (120141 ). Thus, the (two-sided) mass-loss rate 
would be Mj ^ 1.1 x 10 -6 M 0 yr -1 . Since the accre¬ 
tion jyrtehas been estimated to be ^ 5 x 10 -6 M 0 yr -1 
(|Lee et al ] l20li , the mass-loss rate is estimated to be ^ 


associated with the knots (SK2, SK3, SK4, and SK5) in 

X-wind model (~ 30%, Shu et al. 

200 I) 

the jet. The bow wings are faint in the west, causing them 

model (~ 10%.lKonid & Pudritz 

2000 b. 


and the disk-wind 


structure can also be seen along the jet axis in between 
knots SK3 and SK4 at this velocity, tracing the jet itself. 
SiO emission also shows the similar bow shocks curving 
back to the jet axis (see Figure HJd). This feature has been 
seen in the simulations of a pulsed jet, both unmagnetized 


(IStone fc Normanl I19931 [Biro fe Ragal Il994t IVolker et al.l 


1999) and magnetized ([Stone fc Hardeel l2000f) . On the 
other hand, a pulsed wide-angle wind would produce a 
chain of wide b ow shocks that do not curve back to the jet 
axis ([Lee et al.ll2001 h inconsistent with our observations. 

3.3. Proper Motion 

Since the knots in the jet are well traced by the SiO 
emission, proper motion of the jet can be estimated by 
measuring the position shifts of their SiO emission peaks 
with respect to th ose seen at ^ 6 years earlier in our previ¬ 
ous SMA map in ILee et all ( 2008 ). We first convolved our 
previous SMA SiO map to the resolution of our ALMA 
map and then aligned the two maps with the continuum 
peaks. As shown in Figure [5j the jet pattern is similar 
in the two epochs, in agreement with the pattern being 
propagating with the jet velocity. Only the knots that 
are resolved and have well defined peaks are used for the 
measurement, with their emission peak positions marked 
with the solid lines. The innermost pair of knots, SS and 
SN, are not used for the measurement because they are 
spatially unresolved in our ALMA observations, appear¬ 
ing as jetlike structures. Thus, the mean position shift is 
estimated to be ^ 0"36±0"15. Hence, the proper motion 
is estimated to be ~ 0!'06±0"025 per year, giving rise to 
a tangential velocity of ~ 115T50 km s -1 . Since the jet 
is almost in the plane of the sky, the jet velocity can be 
approximated to this tangential velocity. 

3.4. Jet Density and Mass-Loss Rate 

With better estimated jet velocity and better resolv ed 
CO emission than those reported before in ILee et all 
(I 200 I . we can refine the jet density and thus the mass- 
loss rate. The CO emission of the jet in between knot 
NK1 to knot NK5 appears continuous and smooth (see 
Figure Ht), and thus can be used to derive the mean jet 
density. The mean CO intensity there is found to be ^ 


4. DISCUSSION 

4.1. Bow Shock Formation and Sideways Ejection 

As discussed earlier, the knots in the southern jet are 
seen with small bow shocks, with their bow wings curving 
back to the jet axis. Knots SK4 and SK5 are reasonably 
resolved, allowing us to study the formation mechanism of 
the bow shocks. For knot SK5, the CO emission peak ap¬ 
pears clearly upstream of (i.e., closer to the central source 
than) the SiO emission peak (Figure [ 6 ^). For knot SK4, 
the CO emission peak appears only slightly upstream of 
the SiO emission peak by 008 or ^ 15% of the beam size, 
and thus higher resolution observations are needed to con¬ 
firm it. In SiO, the position-velocity (PV) diagram cut 
along the jet axis shows an arclike PV structure for each 
knot, with an apex at the velocity of ^ 9 km s -1 pointing 
away from the source, as shown in Figure [ 6 ) 3 . As we move 
away from the apex, either to higher velocity or to lower 
velocity, the PV structure bends slightly toward the source 
direction, forming the arclike PV structure. The CO emis¬ 
sion shows a similar PV structure slightly upstream of that 
of the SiO emission. Note that the CO emission is missing 
from ~ 8 to 10 km s -1 due to a foreground cloud at those 
velocities. In the northern jet, the knots are not spatially 
resolved and thus no clear bow-like structures can be seen 
in SiO and CO. Higher resolution observations are needed 
to study them. 

Arclike PV structure has bee n seen in pulsed jet sim - 
ulations (see e.g., Figure 16a in IStone fc Normanlll993l) . 
In those simulations, even though the jet source is contin¬ 
uously ejecting material with a constant density, a (sinu¬ 
soidal) variation in the ejection velocity prod uces an inter¬ 
nal working surface (IWS) in the jet dBiro fc Ragalll994 ). 
The IWS can appear as a small arclike knot, with its ve¬ 
locity increasing toward its two edg es fr om the j et axi s 
(see Figure [7] and also Figure 10 in ILee fc Sahail |2004[ ), 
producing the arclike PV structure. As its shocked ma¬ 
terial is ejected sideways into the outflow cavity, it forms 
a bow shock, wit h the bow wings curving back to the jet 
axis ([Biro fc Ragalll99 4i), like those of knots SK4 and SK5 
(see Section lT 2 |) . Therefore, knots SK4 and SK5 likely 
trace the IWSs in the jet. A bigger variation in the jet 
velocity would produce faster and more massive sideways 
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ejection, forming stronger bow shocks ([Raga et al.ll2002l h 
such as SKI, NK1, NB1/2 and SB 1/2, and thus bigger in¬ 
ternal shells extending backward to the central source from 
their bow wings, as see n in the pulsed je t simulations in 
iBiro fe Ragal (| 1994(1 and iLee et all (|2Q0lh . 

Sideways ejection hajybeen claim ed in another jet source 
IRAS 04166+2706 ([Santiago-Garcia et alii 112009 ). In that 
jet, a sawtooth pattern was seen in the PV diagram cut 
along the jet axis in CO and SiO, and it was argued to be 
produced by the sideways ejection of the shocked material 
in the IWSs. That jet is highly inclined with an inclina¬ 
tion angle of ^ 45° to the plane of the sky. Viewing the 
IWSs at high inclination angle produces the sawtooth PV 
pattern, as predict ed in the simulated PV d iagram at high 
inclination angle in lStone fe Normanl ([l993). Here the HH 
212 jet is almost in the plane of the sky with a small incli¬ 
nation angle, thus the sideways ejection forms the arclike 
PV structures in the PV diagram. 

4.2. Wide-Angle Wind Component? 

In addition to a collimated jet, an unseen wide-angle 
wind seems to be needed in driving mo lecular outflows 
in the later sta ge of star formation ( Lee fe Ho 120051: 
lArce et ahll2013f) . Is a wide-angle wind also needed in the 
early phase of star formation as in HH 212? Recently, 
a wide-angle flow has been seen in C 34 S and m ay sug - 
gest a presence of a wide-angle wind ( Codella et al.ll20T3i . 
T his wide- angle flow has also been seen earlier in HCO + 
([Lee et al.ll2014l) and it may trace the cavity walls instead. 
In addition, a wide-angle wind, if exists, would produce 
wide bow-shocks and thus wide internal shells ([Lee et al.l 
l200lf ). This is inconsistent with our observations, which 
show highly curved bow shocks and narrow shells driven by 
the knots. Therefore, in this source the wide-angle wind, 
if exists, is not significant enough to affect the structures 
of the bow shocks and the shells. This is not inc onsis¬ 
tent with current jet l aunching models ([Shu et al.l 120001: 
iKonigl fe Pud ritz 2000), which predict a wide-angle wind 
much more tenuous and thus less significant than the cen¬ 
tral jet. In the later phase, the jet is much tenuous, so 
that the contribu tion of the wide-angle wind could become 
more important ([Lee fc Holl2005l: lArce et aPl2013h . 

4.3. Jet Wiggle 

The jet wiggle can be better seen after combining the 
CO and SiO maps of ALMA with the H 2 map, as shown 
in Figure [51 The northern jet is continuous and the wiggle 
is reasonably well defined. The southern jet appears to 
be continuous only up to knot SK2 at ^ —10", and thus 
the wiggle is not so well defined. Comparing the sinuous 
structures in between knots SKI and SK2 to those in be¬ 
tween knots NK1 and NK2, the wiggle in the southern jet 
could be roughly point-symmetric to that in the northern 
jet. Therefore, we first model the wiggle in the northern 
jet and then apply the model to the southern jet assuming 
a point-symmetric wiggle about the central source. 

A possible mechanism to produce a p oint-symmetric 
wiggle is the jet precession ([Raga et al.ll2009l h In this case, 
the amplitude of the wiggle, A , will increase linearly with 
the distance, i.e., 

(i) 


where 2 measures the distance of the jet from the source 
along the jet axis. Here 0q is the half-opening angle of the 
wiggle, so that tan#o gives the rate of the growth in the 
amplitude with the distance. The wiggle appears to be 
semi-periodical and thus can be assumed to be given by 
the following sinusoidal model 

2 7 T 7 

x(z ) = A(z) sin(—+ M ( 2 ) 

where x measures the displacement of the jet perpendicu¬ 
lar from the jet axis. A is the wavelength of the wiggle and 
4 >0 is the phase angle at the source. Based on the wiggle of 
the northern jet, the best parameters are 60 ~ 0.5°±0.2°, 
A - 576+170 (or 2240+400 AU), and 0 O - 0°±50°. Since 
the jet velocity is ^ 115 km s _1 , the period of the wiggle 
is ~ 93 years. As seen in Figure [8^, the model can roughly 
reproduce the wiggle of both the northern and southern jet 
in the inner part, but it seems to predict a larger wiggle 
than the observed in the outer part with \z\ > 20". 

In order to improve the model, the amplitude of the wig¬ 
gle is assumed to first increase linearly with the distance 
and then remain constant at z > zq, i.e., 


J ztan^o if z < zo, 
1 zq tan 60 if z > zq 


(3) 


As shown in Figure [HJd, the model with zo ~ 10" (4000 
AU) can give a better match up to ^ 34" (i.e., 6 cycles of 
wiggle) in the northern jet and up to ^ 38" in the south¬ 
ern jet, where the SiO and CO emission were detected. At 
z = zo , the amplitude of the wiggle reaches the maximum 
value of A m = zo tan#o ~ 0" L roughly the same as the jet 
radiu s, which is also ~ 071 ([Cabrit et al.l 120071: ILee et all 
120081 ). Note that, since the southern jet does not have con¬ 
tinuous structure beyond ~ —10" from the source, further 
observations are still needed to confirm that the wiggle of 
the jet is really point symmetric about the central source. 
Observations at higher angular resolution are also needed 
to refine the value of zo. 

If the wiggle is really point-symmetric about the central 
source, then one way to produce the wiggle is the jet pre¬ 
cession due a tidal intera ction of a binary com panion on 
the disk of the jet source ([Ter quern et al.lll999l) . However, 
if this is the case, the amplitude of t h e wigg le is expected 
to grow with the distance ([Raga et alJl2009h . inconsistent 
with our observations. 

A nother possibility is t he cur rent-d r iven kink instabil¬ 
ity ([Ceraueira fc de Gouveia Dal Pinoll2QQlh since the jet 
is expected t o be magnetized with a helical magnetic field 
morphology ([Shu et al.l 120001: iKonigl fc Pudritzfl2000l) . In 
this case, the amplitude of the wiggle is expected to grow 
at the beginning, and then become saturated and thus re¬ 
main roughly constant at l arge distance, see , e.g, Figure 
3 in both Cerq ueira fe de Gouveia Dal Pind (2001 ) and 
iMizuno et al.l ([20141 ). In addition, the wiggle is more like 
oscillation and thus could be semi-periodical. Both of 
these are roughly consistent with what we observe here 
in our jet. However, the wiggle is not expected to show 
any symmetry about the central source. This might be OK 
because the jet wiggle here might not be point-symmetric, 
as mentioned earlier. In addition, the detailed jet struc¬ 
ture here already appears to be different between the 
northern jet and southern jet. For the kink instability 
to take place, we have the Kruskal-Shafranov criterion 


A(z) = ztan#o 
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\B p /B ( f ) \ < X/27TA m (e.g., lBatemanlll978t ). With the max¬ 
imum displacement A m ~ 07l and the wiggle wavelength 
A ~ 5"6, we have \B p / < 9. Therefore, this kink 
instability, if in action, could be initiated in the central 
part of the jet, where the magn etic field is dominated by 
the poloidal field ([Pudritz et al.ll2012l) . This poloidal field 
could serve a “backbone” to stabilize the jet ([Ouved et al.l 
l2003l h The toroidal field dominates only near the jet edges 
in order to collimated the jet. 

5. CONCLUSIONS 

We have mapped the HH 212 jet and outflow system 
in the SiO (J=8-7) and CO (J=3-2) lines. Our primary 
conclusions are the following: 

• Although the jet is seen with knots and bow shocks, 
the underlying jet is continuous, especially in the 
northern side. 

• The proper motion of the jet is estimated to be ^ 
115T50 km s -1 . The mass-loss rate in the jet is 
refined to be ^ 1.1 x 10 -6 M© yr -1 , as expected 
for a Class 0 source. 

• Some knots are seen associated with small bow 
shocks, with their bow wings curving back to the 
jet axis. Two of them, e.g., knots SK4 and SK5, are 
reasonably resolved, showing kinematics consistent 
with sideways ejection, likely tracing the internal 
working surfaces formed in a collimated jet due to 
a temporal variation in the jet velocity. Their bow 
shocks are formed by this sideways ejection. 

• The jet has a small semi-periodical wiggle, with a 
period of ~ 93 yrs. The amplitude of the wiggle first 


increases with the distance and then stays roughly 
constant. A jet precession may have difficulty to ex¬ 
plain this kind of wiggle. One possible origin of the 
wiggle could be the kink instability in a magnetized 
jet. 

• Nested (internal) shells are seen in CO connecting 
to the knots and bow shocks, driven by them. The 
sideways ejection ejects material into the outflow 
cavity from the knots, forming the bow shocks and 
the internal shells. 

• The internal shells of the knots are narrow and the 
bow shocks of the knots curve back to the jet axis, 
both suggesting that the wide-angle wind, if exists, 
is insignificant as compared to the jet. 
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Fig. 1.— CO channel maps (Red, in logarithmic scale) in comparison with the H 2 map (Green, adopted from lMcCaughrean et ciD l2002) of 
the HH 212 system. The star marks the central source position. The channel velocity is indicated in the lower left corner. The dashed arcs 
outline the primary beam coverage in our ALMA observations. The synthesized beam in the CO maps has a size of 0"56x0"47 with a PA ~ 
20°, as shown in the lower right corner in (a). 
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Fig. 2.— CO maps (red) integrated over selected velocity ranges that show the shell structures, in comparison to the H 2 map (green). In 
(a), the velocities are from —3.6 to —2.9 km s _1 , from —2 to —1.3 km s _1 , from 1.7 to 2.9 km s —1 , and from 4.9 to 6.1 km s -1 . In (6), the 
velocities are from —7.8 to —5.6 km s _1 , from —3.6 to 0.25 km s _1 , and from 1.6 to 1.8 km s -1 . Nested CO shells are seen connecting to the 
H 2 bow shocks and knots. The dashed lines indicate the connections of the wide-opening shells to the big bow shocks, SB3 and NB3. 
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Fig. 3.— CO channel maps (Red, in linear scale) in comparison with the H 2 map of the HH 212 system in the inner part. The synthesized 
beam in the CO maps has a size of 0"56x0"47, as indicated in the lower right corner in (a). The star marks the central source po sition. The 
channel velocity is indicated in the lower left corner. The contours show the continuum map at 350 GHz adopted from [Lee et al ] usa). 
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Fig. 4. — The jet in SiO, CO, and H 2 . Gray image shows the H 2 map. In ( b ), contours show the SiO map (integrated from —23.03 to 
17.53 km s -1 ). Contour levels are 0.23 x 2.5 n-1 Jy beam -1 km s -1 , where n = 1,2,3... In (c), contours show the high-velocity CO map 
(integrated from —21.65 to —8.95 km s -1 on the blueshifted side and from 9.56 to 18.77 km s -1 on the redshifted side). Contour levels are 
0.20 X 2 n-1 Jy beam -1 km s -1 , where n = 1,2,3... Panels (d) and (e) are the zoom-in of panels ( b ) and (c), respectively. 
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Fig. 5.— Proper motion measurement using the SiO knots in the jet in two different epochs at the ALMA angular resolution of 0056x0047. 
Contour levels are 1.2 x 1.6 n—1 Jy beam -1 km s 1 , where n = 1,2, 3... The star marks the central source position. The line segments indicate 
the SiO peak positions of the knots. 
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Fig. 6.— Structure and kinematics of knots SK4 and SK5 in the south. Gray contours with image are for CO and black contours are for 
SiO. (a) The SiO map is the same as in Figure [4] but zooming into knots SK4 and SK5, highlighting their bowlike structures. The CO map 
shows the total CO emission of knots SK4 and SK5 integrated from 2.2 to 15.4 km s -1 . Contour levels are 0.25 X 2 n-1 Jy beam -1 km s -1 , 
where n = 1,2,3... The dotted line indicates the jet axis in the south. ( b ) Position-velocity (PV) diagrams cut along the jet axis for knots 
SK4 and SK5. The vertical dashed line indicates the systemic velocity of 1.7 km s -1 . The contour levels start at 40 mjy beam -1 with a step 
of 160 mjy beam -1 for CO, and start at 30 mjy beam -1 with a step of 120 mjy beam -1 for SiO. 



Fig. 7. — Velocity structure of the shocked material in an internal working surface (IWS) in the shock frame moving down the jet axis. The 
vectors show the material motion, with the length showing the magnitude. Two layers, CO and SiO, are plotted as an illustration based on 
our observational results. 
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Fig. 8.— The wiggle of the jet and the model (red curves) for the wiggle. Gray image is H 2 , green image is high-velocity CO, and red 
image is SiO, all from Figure [4] (a) The model assumes a constant growth in the wiggle amplitude with the distance from the central source. 
(b) A revised model assumes a constant wiggle amplitude for a distance greater than 10 // from the central source. 















